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We consider polariton condensation in a generalized Dicke model, describing a single-mode cavity 
containing quantum dots, and extend our previous mean-field theory to allow for finite-size fluctua- 
tions. Within the fluctuation-dominated regime the correlation functions differ from their (trivial) 
mean-field values. We argue that the low-energy physics of the model, which determines the pho- 
ton statistics in this fluctuation-dominated crossover regime, is that of the (quantum) anharmonic 
oscillator. The photon statistics at the crossover are different in the high- and low- temperature lim- 
its. When the temperature is high enough for quantum effects to be neglected we recover behavior 
similar to that of a conventional laser. At low enough temperatures, however, we find qualitatively 
different behavior due to quantum effects. 
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I. INTRODUCTION 

Microcavity polaritonsiiSi&i^ are quasiparticles which 
form in wavelength-scale optical cavities containing di- 
electrics. They are mixed modes formed from cavity 
photons and dielectric excitations such as excitons. Since 
they are part photon polaritons are bosons, and are thus 
candidates for Bose condensation. 

Experimental results on pumped microcavities con- 
tinue to be linked to polariton condensation 6 * 7 *^. The 
basic result is a threshold behavior of the luminescence 
intensity from a driven microcavity, while other features 
seen include non-thermal correlation functions for this lu- 
minescence, along with spatially and spectrally localized 
emission. 

The central characteristic of Bose condensation is the 
generation of many-particle coherences which, in polari- 
ton condensation, appear in the electromagnetic field. 
The existence of polariton condensation remains contro- 
versial because polariton condensation is not the only 
phenomenon we can associate with coherent photons in 
microcavities. Most straightforwardly, if the cavity is 
driven into the weak-coupling regime one expects conven- 
tional lasing, and this is thought to be the correct inter- 
pretation of early claims for polariton condensationiiaii. 
While the more recent experiments cannot be straightfor- 
wardly attributed to conventional lasing, more exotic al- 
ternatives to polariton condensation, such as polaritonic 
lasing 12 , remain. 

While there may appear to be many different routes 
to optical coherence in microcavities the relationships 
amongst these routes are not clear. Polariton conden- 
sation, polariton lasing, and conventional lasing are of- 
ten assumed to be fundamentally distinct, but could 
equally well be related phenomena in different parameter 
regimes. This view is supported by recent work showing 
that adding decoherence processes to a mean-field the- 
ory of polariton condensationi^iiii^ leads to a crossover 
from condensation to conventional lasing 1617 . It is also 
suggested by the fundamental connections between equi- 



librium and non-equilibrium phase transitions. The best 
known of these connections is between laser theory and 
the Landau theory of second-order phase transitionai&iS, 
but we note also recent work connecting the critical be- 
havior of parametric oscillators and ferromagnets 20 , and 
a treatment of the ideal Bose gas along the lines of laser 
theorjiSi. 

In this paper we investigate fluctuations close to polari- 
ton condensation, and how they affect the photon statis- 
tics. In general there are different regimes for the domi- 
nance of fluctuations: in a very large system at low den- 
sity of excitation the thermal equilibrium transition is of 
the BEC varietyj 2 -* 2 ^ so that spatial fluctuations are im- 
portant. But since the polariton mass is very small (be- 
cause of the large ratio between the wavelength of light 
and the typical exciton radius and exciton separation), 
at modest densities a BCS-like mean-field regime occurs 
and spatial fluctuations are small. Here the dominant 
fluctuations may be due to the finite size - or generi- 
cally the finite population - of the system. This is very 
often also the case for conventional lasers**^ (for similar 
reasons). 

Here we consider finite-size fluctuations in isolation by 
studying a single-mode model microcavity. We obtain an 
approximate form for the free energy of the model by ne- 
glecting quantum effects. This form can be interpreted as 
the classical probability distribution for the intensity of 
the cavity field, and used to obtain all the static correla- 
tion functions of the cavity photons. It is identical to the 
accepted form for the intensity distribution near the on- 
set of lasing, so we argue that within our approximations 
polariton condensation and lasing are not distinguished 
by the qualitative behavior of the static correlation func- 
tions. However the parameters in the intensity distribu- 
tion are associated with different physics in the two the- 
ories, so there remains room for quantitative distinctions 
between them. 

While neglecting quantum effects in a theory of con- 
densation leads to the same intensity distribution pre- 
dicted by classical laser theory, this approximation fails 
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for a condensate at low temperatures. When the temper- 
ature becomes comparable with or less than the interac- 
tion energy of two photons the finite level spacing affects 
the correlation functions. We shall see that this leads to 
behavior for the correlation functions of a condensate at 
low temperature which is qualitatively different from the 
predictions of standard laser theory. This difference is 
not surprising because standard laser theory is a classi- 
cal approximation, controlled by the photon number at 
threshold. 

Our analysis concerns the thermal equilibrium of a sim- 
plified model of a microcavity. Although in some param- 
eter regimes current experiments are far from thermal 
equilibrium, there are several reasons to study the equi- 
librium behavior. The above-threshold luminescence in 
some recent experiments^ is suggestive of thermal equilib- 
rium, and as microcavities continue to develop—^ exper- 
iments can be expected to reach states closer to thermal 
equilibrium. Furthermore, the behavior close to equilib- 
rium is expected to be similar to that in equilibrium, and 
an understanding of the equilibrium physics provides the 
basis for developing non-equilibrium theories. Finally, 
the qualitative behavior we exhibit here, in particular 
the functional forms of the correlation functions, derives 
from the structure of the effective theory describing the 
collective variables. This structure may be independent 
of whether the theory describes an equilibrium system 
such as a conventional condensate or a non-equilibrium 
system such as a conventional laser. 

The remainder of this paper is organized as follows. 
We begin with an outline of the main results in section 
ITU which compares the classical laser, the classical and 
quantum critical fluctuations of the polariton condensate 
in terms of the anharmonic oscillator. The rest of the 
paper provides a derivation of these results, and a deeper 
quantitative analysis of the model. 

In section IIIII we present the model we consider for 
both lasing and condensation, and give some necessary 
background on the mean-field theory of polariton conden- 
sation. Section Hvl contains the general analysis resulting 
in the free energy in the transition region, including a 
discussion of the regime of applicability of the classical 
approximation. In section^we use these general results 
to analyze the phase diagram and photon statistics of the 
condensate. In section IVT1 we briefly review the standard 
calculations of the intensity distribution in a laser near 
threshold, and compare with our results for the polariton 
condensate. In doing so, we note that conventional laser 
theory and its approximation to a phase transition relies 
on a "large-N" justification that is not usually exposed. 
Section IVIII contains numerical estimates of the size of 
the fluctuation-dominated and quantum regimes in cur- 
rent condensation experiments. In section IVIIII we dis- 
cuss prospects for systems with large quantum regimes, 
and the role of spatial fluctuations. Finally, section IIXI 
summarizes our conclusions. 



II. OUTLINE OF MAIN RESULTS: THE 
ANHARMONIC OSCILLATOR 

We shall find that the classical laser and quan- 
tum/classical condensate lie in the same universality class 
as the anharmonic oscillator 

H = a^(p + j^^(f>(p. (1) 

We choose 4> to be normalized to obey the canonical Bose 
commutation relation [(/>, <j>' i ] = 1. 

The eigenstates of Q are just the number states, \n), 
with energies E(n) = (a — 7)71 + jn 2 an + ^yn 2 . Thus 
the partition function is 

00 

z= ^ e -%»V) j (2) 

with (3 = l/k B T. 

The parameter a is the tuning parameter through the 
transition. At the mean-field level, minimization of the 
exponent in (0) leads to 

"■min =0 for a > 

(3) 



Expanding in terms of fluctuations, Sn = 11 — n min , one 
obtains (on the "condensed" side a < 0) quadratic num- 
ber fluctuations controlled by j(Sn) 2 . We shall choose 
the parameter 7 cx 1/N, with N growing with the sys- 
tem size; the limit N — ► 00 gives the mean- field result. 
We show below that the Dicke model indeed gives rise to 
a partition function of the form of after truncating 
higher order terms in the exponent that are systemati- 
cally smaller in powers of 1/N (where in this case N is 
the number of quantum dots in the cavity). 

The summation in the partition function is of course 
still over discrete quantum states \n), but if the temper- 
ature is high enough the discreteness is irrelevant and 
the sum can be replaced by an integral. We now use 
fields ip = <p/VN rescaled by system size so that ip ap- 
propriately describes the classical electromagnetic field 
intensity. Then one may write the partition function as 

Z w J di,dre- pN[aW2+ ^ NW ' ] (4) 

(remembering that jN = 0(1)). We derive this form ex- 
plicitly for the Dicke model in section lTvl We now remark 
that the action in 10} is consistent with the steady-state 
distribution from a Fokker-Planck equation for diffusion 
in a quartic potential (see section EJ, which is the con- 
ventional approach to laser theory. Of course the "tem- 
perature" there is a fiction that is generated by couplings 
to (Markov) baths representing the outside world. Fur- 
thermore there is no quantum limit for that type of laser 
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theory. In the classical limit the steady state laser and 
the polariton condensate have the same scaling form. 

The anharmonic oscillator thus encapsulates the basic 
results of the two different models when we focus on the 
low energy physics. One signature of the transition is 
the intensity- intensity correlation function g^ 2 \0), which 
is straightforward to calculate for the anhmarmonic os- 
cillator: 



ff (2) (0) 
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which is of course a step function, 
chosen 7 = 1/N = 1 /40 for this 
the qualitative evolution will be 
values of 7. 

While at high temperatures (3> 
classical form, it is very different 
7. Well below the transition the 
to <?( 2 )(0) describes the thermal 
two photon states, so 



For specificity we have 
demonstration, though 
the same for different 

7) g^ (0) resembles the 
at low temperatures, < 
dominant contribution 
excitations of one and 
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Although in this regime the statistics are strongly sub- 
Poissonian, the intensity is also very low, w e~^ a . On 
crossing the transition the intensity increases and (0) 
approaches 1. At the lowest temperature g^ 2 \0) is a 
small constant below the transition, which then approx- 
imately follows the zero-temperature result above the 
transition: g^(0) = l — l/n with n the nearest integer 
to —0/(27) — 1/2. At higher temperatures the quantum 
corrections decrease, and the form of g^ffi) is intermedi- 
ate between the classical result and the low-temperature 
one. 

While the results are straightforwardly exposed in 
terms of the anharmonic oscillator, the mapping of the 
polariton condensate to the parameters a and 7 is more 
complex. The essential details are exposed in the phase 
diagram of Figure [5] 



-a 

FIG. 1: g^(0) as a function of a for the thermal equilib- 
rium of the anharmonic oscillator Q. The top panel shows 
the predictions of the classical approximation, while the bot- 
tom panel is the result of numerically evaluating the partition 
function. 7 = 1/40, and T = 2, 1/2, 1/5, 1/15, and 1/50 for 
the solid, dashed, dotted, dot-dashed, and long-dashed curves, 
respectively. The upper family of lines in the top panel are 
at 4(7/3) = (a/3) 2 , indicating the boundary of the classical 
fluctuation-dominated regime. The lower family are a = T, 
indicating the range of validity of the classical approximation 
in the normal regime. The circles on the lower panel are the 
(discrete) values taken by g^ 2 '(0) at T = 0. 

In Fig^we plot g^ 2 \0) calculated numerically for ||TJ 
as a function of a, at several temperatures. The tuning 
parameter a is proportional to density (at fixed temper- 
ature) for the polariton problem, and generically is the 
pump rate in the laser model. The mean-field ordered 
phase is to the right, where we recover asymptotically 
g^ (0) — > 1 as expected for a classical macroscopic field. 
In the top panel the classical (high-temperature) result 
is shown, with different curves corresponding to differ- 
ent temperatures. As temperature is lowered the result 
approaches more and more closely the mean-field theory, 
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FIG. 2: Fluctuation-dominated regimes in terms of density 
and temperature for the polariton condensation transition in 
the model Q, with E(i) = uj so that all excitons are reso- 
nant with the cavity mode. The solid line is the mean-field 
boundary, and dashed(dotted) lines mark the boundaries of 
the fluctuation-dominated regions for N = 10(40). Curves 
correspond to those in Fig. [3] The shading indicates the re- 
gions in which the static approximation holds for N = 10. 
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The solid line is the mean-field phase boundary be- 
tween the uncondensed (to the left) and condensed 
phases. Here the polariton density is measured per site, 
which explains why the mean-field boundary asymptotes 
to p cx = 1/2 (the boundary of inversion) at large tem- 
peratures. The dashed and dotted lines mark regimes 
on either side of the mean-field transition where we es- 
timate fluctuations to be substantial. These lines cor- 
respond to the upper panel of lines in Figure ^ The 
strong asymmetry of the fluctuation regime about the 
mean-field line arises because we have chosen to expose 
the results with density as a parameter, rather than an 
external tuning parameter (here, naturally, the chemi- 
cal potential). On the condensed side of the transition 
the number fluctuations are large because the chemical 
potential for polaritons becomes nearly clamped. The 
shading marks the estimated regime of validity of the 
classical approximation (akin to the lower panel of lines 
in Figure 0. Over most of the figure the shaded regime 
encloses the fluctuation regime, and the classical approx- 
imation (which guarantees that g^ 2 \0) > 1) holds. For 
small densities and temperatures, however, we see that 
there is a strongly-interacting quantum regime. 

III. MODEL 

Mean-field theories of polariton condensation have now 
been developed from models of a range of systems, in- 
cluding those with propagating photons 22,23 . Here we 
consider the simplest model which leads to a theory of 
polariton condensation, the generalized Dicks 2 ^ model 

ff = ^+£^^ + ^E(^+^r)- (7) 

i * i 

This is the basis of standard laser theory, as well as of 
our earliest mean-field theoryi 3 *!^^ of polariton conden- 
sation. It directly describes a system of localized exci- 
tons, for example in quantum dots, embedded in a three- 
dimensional microcavity. There is a single cavity mode, 
with annihilation operator tp and energy u>, dipole cou- 
pled to a set of N quantum dots or atoms of the gain 
medium. The state Sf — +(— )1 corresponds to the pres- 
ence (absence) of an exciton on site i, or to an atom 
in the upper (lower) of the states coupled by the lasing 
transition. 

We have included a factor of \/\/~N explicitly in the 
light-matter coupling in JIJ), so that for a cavity mode 
coupled to N dots g is related to the observed Rabi split- 
ting and not to the single-dot Rabi splitting. This scaling 
is formally necessary because we will be concerned with 
the properties of the model (0 for large N, so need the 
model to be well behaved in the limit ./V — > oo. But we 
stress that with our convention changing the number of 
dots in a cavity of fixed volume corresponds to changing 
both N and g. 

To construct a theory of Bose condensation of polari- 
tons from Q one studies its thermodynamics, fixing the 



total number of excitons and photons 

^Vcx = + \ E i S i + 1) ■ ( 8 ) 

i 

Thus we consider the free energy, which can be written 
as the functional integral 

JL=f = \nJvi/ J Vje- NS °«, (9) 

where 

S cS = / dTi/j(d T + £j c )4> - lndet P. (10) 
Jo 

Here if> is related to the real electromagnetic field ampli- 
tude ipo by ipo — ipVN, while lndet P corresponds to the 
free energy of the quantum dots in the field ip(r). 

Eq. (0 is the free energy in the grand-canonical en- 
semble. The constraint on N ex is dealt with on aver- 
age, by introducing a chemical potential /i which appears 
in (|10|l as a shift of the photon and exciton energies: 
uj — ► Hi = (w — fi) and E(i) — > E(i) — (E(i) — //,) = 2ii. 
In the limit N — ► oo the relative fluctuations of N ex 
which occur in the grand-canonical ensemble are negli- 
gible, and the grand-canonical and canonical ensembles 
are expected to be equivalent. However, we are now con- 
sidering fluctuations in a finite system, for which the pre- 
dictions of the canonical and grand-canonical ensembles 
could differ. The correct ensemble then depends on how 
the polariton number is constrained experimentally. We 
consider an idealized limiting case in which this is done 
by coupling to a reservoir, so that the grand-canonical 
ensemble is appropriate. 

We can develop an asymptotic expansion of the free 
energy for large N by expanding S c g around its static 
saddle points ip( T ) — 4>- The saddle-point equation gives 
the mean-field estimate of the phase boundary between 
the normal state (ip = 0) and the condensed state (ip 
finite). Evaluating @ on the stable saddle point leads 
to the mean-field estimate of the free energy; this is the 
only extensive contribution to the free energy, so that the 
mean-field theory gives the exact asymptotic form for the 
free energy in the thermodynamic limit N — oo. 

Polariton condensation in the model can be viewed 
as a generalization of the ferroelectric transition discov- 
ered by Hepp and Lieb—^. There has been recent inter- 
est in the physics associated with that ferroelectric tran- 
sition in a finite model at zero temperatur e) 28 ! 29 . This 
differs from the present problem in that there is no con- 
straint on the polariton number, i.e. \x = 0, and as a con- 
sequence the rotating-wave approximation implicit in Q 
cannot be made. This leads to qualitative differences^ 
in the dynamics of the unconstrained models at T = 0. 
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IV. GENERAL FORM OF THE FREE ENERGY 
NEAR THE TRANSITION 

Away from the mean-field phase transition the saddle- 
point expansion provides a systematic approximation 
scheme for the free energy of the finite system. How- 
ever, this scheme fails in the vicinity of the mean-field 
transition due to the diverging occupations of the soft 
fluctuations. In this section, we shall calculate the gen- 
eral form of the free energy of a large but finite system 
near the mean-field transition. To simplify the notation 
we shall take all the excitons to have the same energy, 
E(i) — Eg, but the resulting forms may be straightfor- 
wardly generalized to allow for a distribution of exciton 
energies. 

We consider the action obtained by expanding (|1(J[I to 
fourth order in 



S cS = So + ^X^M^M 



4 £ 



[iuj + E+)(iw + E_) 
(iw + 2e) 

^OJi 1pcU2 T^4 VuJiLJ2UJ3UJ4 



(11) 



So is the action of the two-level systems in the absence 
of photons, while the remaining part describes photons 
in the medium of the two-level systems. There are res- 
onances at E±, which are the polariton energies relative 
to the chemical potential 
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E± = -(w c + 2s± y/(u c - 2s) 2 + 4 ff 2 tanh (J3e)), (12) 

while P1234 is the photon-photon interaction mediated by 
the two-level systems, 



V1234 = y^{(^/ - g) 1 [i{u>f + u 3 ) + £]~ 



x [i(uf + lo 3 - u 2 ) ~ e\ 1 
x [i(uif + UJ4 + uj 3 - 0J2) + e]^ 1 }. (13) 

lui, etc. are bosonic Matsubara frequencies, while ujf — 
(n + 3/4)2irT is a fermionic frequency, shifted to take 
account of the two-level constraint 23,30 . 



A. Static free energy 

As we approach the transition from the normal side 
one of the polariton energies E± is approaching zero, and 
perturbation theory fails. Over most of this fluctuation 
regime, however, u>x — 2nT is large compared with the 
energy of the soft mode. The dominant contribution to 
the free energy for a large system close to the transition 
then comes from the static (classical) paths. Retaining 



only these dominant contributions gives an approxima- 
tion for © near the transition, 



/ » In / d^dre- N{ ^ w2 ~ lacosh ^ i2+! > 2W2) (U) 



= /o+ln / d^*e-"<°M" + WI 4+c M 6+ ->. (15) 



Here a = (3E+E- / '(2e) is the static part of the Gaussian 
kernel in (|llf> . measuring the distance to the transition, 
b is the static part of the interaction term, and c,d, . . . 
are higher-order interaction strengths which do not de- 
pend on the system size N. fo — A^lncosh/Je is the free 
energy of the two-level systems in zero field. The explicit 
form for b, obtained either from i|13|) or by expanding the 
exponent in (|14fl . is 



b = g 4 p 



tanh(/3e) — (3e sech 2 (fie) 
8£* ' 



(16) 



Reduction to an oscillator 



At low temperatures the non-perturbative regime may 
lie outside the regime of validity of the static approxima- 
tion, so that time-dependent paths must be considered. 
The divergences of perturbation theory are still at small 
frequencies, so for a large system close to the transition 
we can replace the action with its low frequency form. 
Considering for definiteness the region near the transition 
where E_ vanishes, we have as small parameters, 

while generically E + and e are finite. The Gaussian term 
in Ijlll) can then be straightforwardly approximated as 



0E. 



2s 



(17) 



with corrections proportional to the ratios of small to 
finite parameters. Approximating the interaction term 
is more involved, because the result of the summation 
(|13|l takes different forms depending on how many of 
the external frequencies coincide. This problem can 
be avoided by restricting our attention to the low- 
temperature regime T -C £, in which all the forms lead 
to the same low-frequency approximation 

8e 3 



E 



Ipl 1p 2 ^3 ^4 



(18) 



^1+^2—^3+^4 



The corrections to (I18f) are again small in terms of the 
ratios of small to large parameters, e.g. oj/s and T/e. 
Since the coefficient in <|18fl is the low-temperature limit 
of 6, it is convenient to replace it with 6, leading to a 
low-energy approximation to the original theory which is 
valid at the quantum level for low temperatures, and at 
the classical level at higher temperatures. 

H17|) and (|18fl become, after rescaling the fields ip — > 

ip.f F 2e N , the action for a quantum anharmonic oscillator 



JTJ . The oscillator frequency a = £L_ , and the interaction 

2 4 



strength is ^ (j^-J 



2iNE\ 



at low temperatures. 
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C. Validity of the static approximation 

The static approximation leading to Ijl4(l holds because 
the energy of fluctuations is much less than temperature. 
On the normal side of the transition this only occurs in 
a region close to the transition. However, the behavior 
of the excitation spectrum on the condensed side is dif- 
ferent. This spectrumiiiiii^ comprises a mode which is 
always at zero frequency relative to /i, and two modes 
which at the transition are at positive and negative of 
the non-vanishing polariton energy. Thus there appears 
to be no mode which becomes soft at the transition, and 
so the classical approximation does not appear to be con- 
trolled by the distance from the transition. 

This puzzle is resolved by inspecting the anharmonic 
oscillator J5J), which we introduced earlier in section ITU 
On the normal side of the transition the nonlinear terms 
in J5J are irrelevant, and the parameter controlling the 
classical approximation is (3a. On the condensed side the 
exponent in @ can be rewritten in terms of the number 
fluctuations Sn = n—n min = n—(—a/(2j)) as — (3^f{8n) 2 . 
Thus the parameter controlling the static approximation 
on the condensed side is f3j. Returning to the polariton 
model close to the transition where EL vanishes, we see 
that the classical approximation holds for the low-energy 
fluctuations when 



b_ 

N 



2s 



<1, 



(19) 



provided the temperature is small compared with the 
Rabi splitting and e. At higher temperatures we expect 
an inequality of the same general form, TN/g > 1, but 
with numerical differences due to the renormalization of 
the effective interaction by the occupation of the high- 
energy polariton and the frequency dependence of the 
interaction. 

Our previous computationsi&iaii^ of the fluctuation 
spectrum were done at a Gaussian level, i.e. approximat- 
ing the partition function of the fluctuations with that 
of a harmonic oscillator. This approximation predicts 
that the characteristic frequency of fluctuations about 
the condensate is zero because there is no linear term 
in Sn when the exponent of J2J) is expanded around a 
finite n. To obtain a finite level spacing one must go 
beyond the Gaussian approximation and include interac- 
tions, which generate a finite level spacing ~ l/N. This 
may be contrasted with the normal state, where the har- 
monic oscillator part of (J2J is enough to generate a finite 
level spacing a. 



CORRELATION FUNCTIONS AND PHASE 
DIAGRAM 



The correlation functions of the cavity field can be ob- 
tained by differentiating l|14|) with respect to uj c or 115|) 
with respect to a. Note that the field is classical: the in- 



tegrand of l|14|) can be interpreted as a probability distri- 
bution for the intensity of the cavity field. Non-classical 
fields are associated with time-dependent paths, which 
give complex integrands in the path integral J^J that can- 
not be interpreted as classical probabilities. 

The integral l|14|l is only tractable numerically. How- 
ever, the interactions c,d, ... do not affect the asymptotic 
behavior of the correlation functions as N — > oo if a > 
and b > 0, because the quartic nonlinearity restricts fluc- 
tuations of the field to \ip\ < N^ 1 ^^ 1 ^ . Therefore we 
can obtain the asymptotic forms of the correlation func- 
tions on the normal side of the mean-field phase bound- 
ary from (|15fl with c, d, . . . = 0. The corresponding free 
energy is, discarding additive constants, 

/ = In / d^dre~ N{aW2+bw4) (20) 

a I N ' 



In- 



e^erfc fJf 



(21) 



On the condensed side where a < we may still use 
the approximation H20(l so long as the order parameter 
\ip | is small. It is a weaker approximation than on the 
normal side, because the minimum of the exponent in 
<P~5l> occurs for |^| 2 ~ 1, so that truncation to a quar- 
tic theory produces errors in the leading asymptotics of 
the correlation functions. However, these errors are pro- 
portional to the order parameter, and so are numerically 
small close enough to the transition, even if they are not 
asymptotically small in N. 

We now present explicit results for the behavior of the 
cavity field in the model © when N is finite but large. 
For orientation, the mean-field phase boundaries E± = 
are shown as the solid lines in Fig. [21 p e x is the number of 
polaritons per site, (N ex /N) with iV ox given by JSJl; note 
that this differs from the definition of p ex used in Refs. 
ITMT4IIT5I by a shift of 0.5. Temperature is expressed in 
units of g, which is one-half of the collective Rabi splitting 
at resonance(see Eg. 11211 . 

For completeness we begin by considering the region 
far from the transition, where the correlation functions 
can be obtained using the saddle-point expansion. On 
the normal side of the mean-field transition there is no 
saddle-point contribution to the photon density (ip^ip), so 
the first non- vanishing contribution is at order l/N. This 
term is just the expectation value of the photon density 
from the Gaussian part of i|ll|) . i.e. the photon density in 
a population of non-interacting polaritons with energies 
E±: 



^ = (E+ - 2e)n B (E + ) + (2£ - E_)n B (E_) 



N(E^ 



E_) 



n B (x) 



a/3x _ ^ ' 



(22) 
(23) 



Since at this order in N we have non-interacting particles 
the many-photon correlation functions are related to the 
photon number by Wick's theorem. In particular, for the 
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static part of the correlation function measured by Deng 
etal. wehave<?( 2 )(0) = (^VW'VW'W = 2+0(l/N). 

On the condensed side of the mean-field transition the 
leading-order contributions to the correlation functions 
come from the saddle point. Thus (ip^ip) ~ 1, and we 
expect 5 (2) (0) = 1 + 0(1/N). Calculations of the higher- 
order terms are complicated due to the presence of the 
zero mode, and so we shall not pursue them here. 

In the region close to the mean-field transition the 
saddle-point expansion must break down, to allow the 
correlation functions to smoothly interpolate between 
their forms in the two states. We can estimate the bound- 
aries of this crossover region by equating the magnitudes 
of successive terms in the large-N expansion of the pho- 
ton density, (i/j^i/j). For the approximate form 1)21(1 we 
see that the crossover regime obeys Na 2 < 46, or 

p N % +E 2 (tanhGSe) - /3esech 2 (/3e)) > 1. (24) 

Note that the numerical factor in this expression is not 
meaningful, since it depends on the precise definition of 
the boundaries of the crossover regime. 



LO — 



CD 



CM — 



O 




5 -0.3 
(u-E g )/(2g) 



-0.1 



FIG. 3: Fluctuation-dominated regimes for the polariton con- 
densation transition with A = 0. Solid lines are the mean- 
field boundaries. Dashed (dotted) lines mark the boundaries 
of the fluctuation-dominated regions for N = 10(40). 

In Fig. we show the phase diagram of the system 
as a function of fi and T, for Ag = to — E g = 0. The 
solid lines are the phase boundary of the infinite system, 
E + E- = 0, while the remaining lines are the boundaries 
of the crossover regions 1(24(1 for systems of 10 and 40 
quantum dots. One general feature which can be seen 
in this figure is that the size of the fluctuation regions 



scales as 1/y/N, as can be deduced from l(2T]l. The par- 
ticular shapes of the fluctuation regions regions come 
from the interplay between the thermal occupations and 
the strength of the interaction, which is temperature and 
density dependent. In particular, the present model has 
the unusual feature that the fluctuation regions vanish 
in the high-temperature limit, where the limits of the 
inequality l(24|) are approximately 



4e = g 2 f3-2Ag± 



Wg 2 



(25) 



Physically, this occurs because for large T the number 
of fluctuations at fixed chemical potential behaves as 
nj uc = ub{E + ) + ub(E^) ~ T while the interaction 
strength b vanishes like 1/T 4 . Thus the mean interaction 
energy (n\ nc b) vanishes in the high temperature limit. 
The occupation of the fluctuations diverges only as T 
because they are confined to a finite energy range, so 
that increasing temperature does not increase the num- 
ber of relevant fluctuation modes, as it would in the more 
familiar case of a semi- infinite band of states. The inter- 
action vanishes because at high temperatures the field 
does not affect the occupation of the two-level systems: 
the free energy of a two-level system, which appears in 
the exponent of ((14(1 . becomes independent of field as 
T — > oo. 

Fig. |2l shows the same phase diagram in terms of den- 
sity and temperature, obtained by relating chemical po- 
tential to density using the mean-field results. For the 
normal state we use 



p cx = Poxcitons = (1 - tanh(/3e))/2, 
while for the condensate the quartic theory gives 

,2 ~ 



Pcxcitons 



(26) 



(27) 



These relations introduce qualitative differences between 
the fluctuation-dominated regions in this figure and those 
in Fig.|5J In particular, the fluctuation-dominated region 
becomes much larger on the condensed side, because the 
chemical potential is only weakly dependent on density 
in the condensate. 

The shading on Fig. [21 indicates the regimes of validity 
of the classical approximation, combining the analysis of 
sections llVBl and iFvTI with Eqs.|23and[23 On the nor- 
mal side of the mean-field transition we shade the region 
[3E- < 1. On the condensed side the analogous inequal- 
ity is ((19(1 . However, the numerical factors in this result 
are only correct close to the phase boundary. Since this 
only occurs at low densities we have plotted the inequal- 
ity T/g > 1/(4N), corresponding to the low-density limit 
of . 

In the upper panel of Fig. 0] we plot the prediction of 
((21(1 for the number of photons in the the cavity 



M 

da ' 



(28) 
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as a function of the deviation of /i from its mean-field 
critical value /i c . We take A = 0, and plot curves for 
T/g = 0.25, 0.75 and 1.5, and for N = 10 and 40. 
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FIG. 4: Photon number N{ip'ip) (top panel) and density 
{ifr ip) (bottom panel) as functions of (/i — fj, c )/(2g) for A = 0. 
T/g = 0.25 (solid lines), 0.75 (dashed lines), and 1.5 (dotted 
lines). N = 10 (lower line of each pair in the top panel, upper 
in the bottom panel) and N = 40 (upper line in top panel, 
lower in bottom panel). Shading marks the condensed region 
of the mean-field theory. 

In general each of the N two-level systems makes a con- 
tribution of order 1/N to the number of cavity photons, 
due to the scaling of the coupling constant. Well below 
the transition we could neglect the interactions between 
the photons and approximate the integrand of (|15f) by a 
Gaussian. This gives Nty'i/j) ~ 1/a. In the Gaussian 
regime each two-level system contributes incoherently to 
the cavity field, so the total photon number is of order 1; 
such scaling is demonstrated by the collapse of the curves 
corresponding to different system sizes in the left side of 
the top panel of Fig. 0] 

As we increase the chemical potential through the 
mean-field transition the occupancy of the cavity field in- 
creases, and the interactions begin to generate coherence 
between the contributions of the different two-level sys- 
tems. Far in the condensed state this coherence is com- 
plete: all the ~ N two-level systems contribute coher- 
ently to the cavity field, giving a photon number which 
scales as N. This scaling can be seen on the right of the 
lower panel of Fig. 0] which shows the photon number 
per two-level system. The order of N contribution to the 
photon number comes from the saddle point of (|20ll , and 
is —a/ (2b). This saddle-point contribution is the expo- 
nential in (|21[1 : it survives in the condensed state, but is 
canceled by the asymptotic expansion of the error func- 



tion well into the normal state. 

In the region near the transition neither the Gaus- 
sian nor saddle-point approximations are appropriate, 
and the full form l|21|l must be used. There is partial 
coherence amongst the two-level systems, leading to a 
photon number which scales as Explicitly, we find 

N(il)'il>) = y/ N/ (tt6) for the photon number at the tran- 
sition. The scaling can be shown from a more general 
argument by noting that the a dependent part of 120(1 is 
a function of y/Na and b, so at the transition 



N(^ijj) 



df(VNa, b) 



da 



a=0 

Ng(0,b). 



(29) 



Fig-EJillustrates the results obtained from (|21|l for the 
static intensity-intensity correlation function 



9 {2, (r = 0) = 1 



Oa- 



2 ' 



(30) 



Well below the transition g^ 2 '(0) approaches the value of 
2 associated with a thermal mixture of non-interacting 
photons, while well above it £P 2 ^(0) approaches 1, corre- 
sponding to the bunched photons of e.g. a coherent state 
or large-amplitude number state. The crossover occurs 
over a range of chemical potentials which scales as 1/ y/N. 



VI. LASER FLUCTUATIONS 

In this section, we briefly review the standard theory 
of fluctuations in a single-mode laser close to threshold, 
following the work of Riskeni&. We shall see that the 
photon statistics predicted by this theory are similar to 
those predicted by the theory of polariton condensation. 
This similarity holds in spite of the fact that the theory of 
condensation describes a system in thermal equilibrium 
while laser theory describes one which is not in thermal 
equilibrium. 

The canonical laser theories are based on the dynam- 
ics of models such as JJJ in the presence of dissipation 
processes. Such dissipation processes can be formally 
generated by coupling each dot, and the photon mode, 
to an oscillator bath. Thus one considers Hamiltonians 
of the general form 



H T = H + H rc 



H s -i 



(31) 



where H = H ph + H s + H s ^ ph is given by Eq. 0, H rcs 
describes the harmonic oscillator baths, and the remain- 
ing two terms couple the system and baths. The stan- 
dard model is to couple the cavity mode to its bath with 
Hamiltonian 



#ph-rcs = X! 9p(^d p + 4V0> 



(32) 
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FIG. 5: Intensity-intensity correlation function (/ 2 ^(0) of 
the cavity photons as a function of chemical potential {(i — 
Hc)/(2g) for A = 0. T/g = 0.25 (solid lines), 0.75 (dashed 
lines), and 1.5 (dotted lines). N — 10 (lower line of each 
pair in the normal region) and N = 40 (upper line). Shading 
marks the condensed region of the mean-field theory. Circles 
correspond to the lines marking the fluctuation-dominated 
regimes of the normal state in Fig. |5] The vertical line corre- 
sponds to the limit of the shading in the normal state on Fig. 
for T/g = 0.25; the analogous lines for T/g = 0.75 and 1.5 
are off the scale. 



and the dots to their baths with Hamiltonian 
H s -rcs = J]] 9pi(dliS i + S^dpi) + g\ p d} pi dpiS z . (33) 

i p 

To develop dynamical equations for the system vari- 
ables alone one takes the Liouville equation for the evo- 
lution of the density matrix and averages over the reser- 
voir variables. The interaction between the system and 
reservoir variables is treated using second-order pertur- 
bation theory, leading to an integro-differential equation 
for the time evolution of the reduced density matrix p. 
This equation is then approximated by a first-order dif- 
ferential equation, implicitly neglecting memory effects, 
formally by having Markovian baths whose spectra are 
completely white. 

Finally, one integrates the atomic degrees of freedom 
out of the density matrix for the system, assuming that 
the atomic system is very strongly damped. This leads to 
an equation of motion for the P representation of the field 
density matrix, W(ip, ip* ,t). Neglecting terms associated 
with the quantum nature of the field 19 this equation be- 
comes the Fokker-Planck equation describing diffusion in 



a quartic potential: 

FfW 8 r) 

d 2 W 

dtp dip* 

The coefficients in this form are the strength of the non- 
linearity ^(denoted in Ref. IT^l. the linear gain or loss 
£ci, and a diffusion constant q. 

To compare with the theory of fluctuations near polari- 
ton condensation we make the same rescaling of the field 
as we did there, tp — > ipvN. In terms of this rescaled field 
the steady-state probability distribution obtained from 
()34|l is proportional to the integrand of (|20|l . The param- 
eters are related to the laser parameters as a = —(d/ (2q) 
and b = £N/(4q). The scaling of these terms follows 
from the microscopic expressions for £, q and d given in 
Eq. (2.17, 2.27) of Ref. [3 Substituting g -> g/VN 
so that g is as defined in Eq. [7| and noting that the to- 
tal inversion a scales as N, we find q ~ N°, ( ~ N^ 1 , 
(d ~ N°. Therefore, as in the polariton condensate, the 
parameters a and b are independent of N. 

Thus in the classical approximation the photon statis- 
tics of a laser and a condensate differ only due to the 
dependence of a and b on the microscopic parameters 
of the system. This is similar to the usual universality 
which occurs close to a second-order phase transition, 
where the singular parts of the observables are given by 
universal scaling functions, with non- universal relations 
between the scaling parameters and physical parameters 
such as temperature. In the finite systems considered 
here the observables are not singular, but for large N 
the correlation functions in the two problems are identi- 
cal functions of a parameter describing the distance from 
the transition. This arises because the collective behavior 
of the electromagnetic field in both theories is described 
by a quartic "free-energy" , and the system size N enters 
these free energies in the same way. 

Within our approximations, then, the only scope for 
making a distinction between polariton condensation and 
lasing arises from the forms of a and b in the two sys- 
tems. These forms are model-dependent, and expected 
to change in more realistic models of either system. How- 
ever, we note that there are differences in the physics of 
these forms in the laser and the condensate, which could 
be expected to be robust. For example, in the laser the 
noise strength q is the spontaneous emission rate into 
the lasing mode, while in the condensate the analogous 
noise strength is temperature. It appears that the laser 
amplifies spontaneous emission noise, while the conden- 
sate amplifies thermal noise. The parameters a and b 
also contain different physics: in the laser a is associated 
with gain or loss and b with gain depiction, while in the 
condensate a is a single-particle energy and b an inter- 
action strength. Thus the size of the threshold regime 
in absolute units is controlled by different physics in the 
two theories. 
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The analogy between the laser threshold and a second- 
order phase transition is well-known, and previous au- 
thors have pointed out that this analogy extends beyond 
the thermodynamic limiliiSi^ 1 : in a finite system the sin- 
gularities associated with a phase transition are rounded 
in the same way that the sharp singularities of the las- 
ing transition in a rate-equation treatment are rounded 
by noise. We see that for the laser model (|31fl the pa- 
rameter controlling this rounding is N, the system size, 
exactly as in the usual thermodynamic case. This ap- 
pears to differ from the picture obtained from classical 
stochastic models 31 , in which the corresponding param- 
eter is the fraction of the spontaneous emission directed 
into the lasing mode. 



VII. QUANTITATIVE ESTIMATES FOR 
CURRENT EXPERIMENTS 

We can use our model to obtain some indication of the 
scale of the fluctuation-dominated region in current ex- 
perimental systems. Such systems are generally planar 
microcavities, rather than the three-dimensional cavity 
of Q. However the luminescence which may be evi- 
dence for condensation is localized, suggesting that the 
condensate itself is localized, presumably due to some 
combination of kinetics, self-trapping, and disorder. In 
a condensate localized on scale L we expect there to be 
an energy cost ~ (L 2 m 2 ) -1 associated with spatial vari- 
ations, where m p ~ 10~ 5 7n e is the in-plane mass of the 
polaritons. If this energy cost is large compared with 
temperature we expect spatial fluctuations to be unim- 
portant, and the present theory to apply. 

Let us consider in particular the data discussed by 
Weihs, Deng, Snoke, and Yamamoto in Ref. HI These 
experiments are done on a GaAs planar microcavity with 
twelve quantum wells and a Rabi splitting of 14.9meV, 
so g w 86K. Their condensate is typically 5/mi in diam- 
eter, giving a temperature scale for spatial fluctuations 
of 60K. We take the Mott density, 10 12 cm -2 , to be an 
upper bound on the density of available exciton sites, so 
the condensate of diameter 5/um corresponds to N ~ 10 5 . 

The temperature in these systems is extremely low 
relative to the Rabi splitting: T — 5K corresponds to 
T/g = 0.06. At these low temperatures the critical den- 
sity is not associated with the single-mode phase tran- 
sition studied here but with spatial fluctuations, whose 
effect on the mean-field theory is analyzed in Refs.l22land 
H3 Thus our theory cannot describe the transition that 
is being crossed in the experiments. However, once the 
condensate has formed we expect it to apply, since the 
condensate which forms is apparently localized, and the 
spatial degrees of freedom are frozen out. We estimate 
that the achieved condensate densities of 10 10 cm -2 cor- 
respond to p cx ps 10 10 /?iMott = 10~ 2 . For this density the 
transition temperature of our theory is T/g = 0.2 w 20K. 
For TV = 10 5 the fluctuation-dominated region extends 
over around 20mK on the normal side of the transition 



and IK on the condensed side. The experimental temper- 
ature is far outside this region, so the present theory sug- 
gests that the fluctuations of the condensate are negligi- 
ble and the light from the condensate area is almost per- 
fectly coherent. If raising the temperature did not change 
the localization then this would allow the sharp crossover 
predicted here to be observed. This complication would 
be absent for systems with external in-plane confinement 
of the polaritons, such as pillar microcavities^ 3 .. 

For a system in thermal equilibrium, we have argued 
that the quantum regime occurs when the temperature is 
smaller than the characteristic level spacing produced by 
the photon-photon interaction. For the particular model 
used here this is when TN ~ g. We estimate that this 
gives a scale of lmK for the systems of Ref. 113, which is 
far below the achieved temperatures. Furthermore, these 
are open systems coupled to baths. In such a system, 
the effect of dissipation is loosely to broaden the level 
spectrum through the decay process. To have quantum 
effects that dominate for the longest times one then re- 
quires linewidths that are narrow in comparison to the 
level spacing. This is not achieved, since a temperature 
scale of lmK corresponds to a linewidth on the order of 
10//eV, far below the linewidth associated with the decay 
of the cavity mode. 



VIII. DISCUSSION 

Recently, three groups^ 3 ^* 3 ^ have achieved strong- 
coupling of single quantum dots in microcavities, and are 
thus beginning to approach the regime, already achieved 
in atom cavity optics 3 ^, of a single two-level system 
strongly coupled to a cavity mode. Clearly the nonlin- 
ear regime 38 of such a system is dominated by quantum 
effects. However, we note that our model suggests that 
this single-atom limit is not the only way to see strong 
quantum effects. For a cavity of fixed volume one should 
replace g/vN with g in ffl. The effective photon-photon 
interaction then scales as ./V(whereas it scales as iV -1 at 
fixed dot density), so that quantum effects survive to 
higher temperatures when the number of dots increases 
in a cavity of fixed volume: a situation opposite to the 
normal thermodynamic limit we consider elsewhere, in 
which quantum effects move to lower temperatures as 
the system size increases. In this context it is interesting 
to note the experiments of Ref. 39, in which the quan- 
tum statistical properties of light emitted from a weakly 
driven atom-cavity system are measured. The observed 
g( 2 )(0) was below the classical limit of 1, and furthermore 
its value was apparently independent of the number of 
atoms in the cavity. This would naturally arise if the 
noise increased as N, as it does in standard laser theory, 
canceling the increase in the level spacing. 

The behavior of a finite system in which spatial fluc- 
tuations are allowed could be analyzed using renormal- 
ization group arguments^. We expect it to depend on 
the interplay between system size and an effective co- 
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ordination number or interaction range: a small system 
with a large co-ordination (long-range interactions) will 
be dominated by finite-size fluctuations over a larger re- 
gion of the phase diagram than it is dominated by spatial 
fluctuations, whereas the reverse will occur for a large 
system with a small co-ordination (short-range interac- 
tions). In the former case we can view the finite-size 
fluctuations as corrections to mean-field theory, as here. 
This cannot be true in the latter case, however, since the 
theory of the finite system should then involve critical 
exponents different from those of mean-field theory. 

The phase diagram of the infinite two-dimensional sys- 
tem, allowing for spatial fluctuations, has been studied 
in Ref. The deviations between the phase bound- 
ary there and that of mean-field theory indicate the re- 
gions in which spatial fluctuations of an infinite system 
are significant. Since the photon mass is very small, ef- 
fectively providing a long-range interaction between exci- 
tons, these deviations are only significant at very low den- 
sities, where even a long-range interaction gives a small 
co-ordination number. Hence, except at low densities, we 
expect even a relatively large system to have mean-field 
like finite-size behavior. 



IX. CONCLUSIONS 

In this paper we have analyzed the behavior of a model 
of a finite polariton condensate close to the mean-field 
transition, and shown that it is formally similar to the 
theory describing a laser close to the laser threshold. This 
similarity fails when the polariton condensate reaches a 
low-temperature quantum regime, since the laser theory 



is classical - notwithstanding that the noise in the laser 
theory might ultimately have its origins in quantum ef- 
fects. 

In the classical regime we find that the intensity distri- 
bution for the photons in a finite polariton condensate is 
of the same form as that obtained from conventional laser 
theory. Thus the photon statistics are not expected to 
reveal any fundamental difference between a condensate 
and a laser, but of course the parametrization is differ- 
ent and the scales have different meanings. In particular, 
for a condensate in thermal equilibrium the sharpness of 
the crossover is controlled by temperature, whereas in a 
laser it is controlled by the noise introduced by coupling 
to external baths. 

Quantum effects dominate in general when the level 
spacing exceeds the noise strength. We have seen that 
for the polariton condensate in equilibrium the relevant 
level spacing is the photon-photon interaction, and the 
relevant noise strength is temperature. In the conven- 
tional thermodynamic limit of our model the tempera- 
ture scale for quantum effects decreases with system size, 
but the opposite occurs when increasing the dot number 
with fixed cavity volume. Hence it is not always the 
smallest systems which have the most strongly quantum- 
mechanical collective behavior. 
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